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^ ; Abstract 

. We apply ideas from renormalization theory to models of cluster formation in nucleation and growth 

o ■ 

O ' processes. We study a simple case of the Becker-D5ring system of equations and show how a novel coarse- 



graining procedure applied to the cluster aggregation space affects the coagulation and fragmentation rate 
coefficients. A dynamical renormalization structure is found to underlie the Becker-Doring equations, nine 
archetypal systems are identified, and their behaviour is analysed in detail. These architypal systems divide 
into three distinct groups: coagulation-dominated systems, fragmentation-dominated systems and those 
systems where the two processes are balanced. The dynamical behaviour obtained for these is found to 
be in agreement with certain fine-grained solutions previously obtained by asymptotic methods. This work 
opens the way for the application of renormalization ideas to a wide range of non-equilibrium physicochemical 



! processes, some of which we have previously modelled on the basis of the Becker-Doring equations. 

o 
o 



I Introduction 

In this paper, we study tlie Becker-Doring cluster kinetic equations familiar from classical nucleation 
theory in which the monomer concentration (ci) is held constant 

Here Cr represents the concentration of clusters containing r-monomers, the dot implies a time- 
derivative, and Jr is the flux from clusters of size r to those of size r + 1. There are certain 
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Figure 1: The effect of our coarse- graining dynamical renormalization on the two parameters 6,p in a 
Becker-Doring model with rate coefficients which vary as a power p of the cluster size r, and with ratio of 
aggregation to fragmentation rates 9. The dots show the fixed points of the mapping. 

mathematical properties of the Becker-Doring system not immediately apparent from the equations 
but crucial to its wide-ranging physical applicability. Firstly, the partition function, Qr, satisfies 
ttrQr = K+iQr+i together with Qi = 1 and formally yields the equilibrium solution c^'^ = QrC^. The 
relevance of this solution depends on the behaviour of Qr in the limit r —>■ oo; further analysis of 
this is given in Section V, where specific examples are analysed in detail. The function 



V{{Cr})=J2cr(\0^ 



Cr. 
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(2) 



r=l \ \Qr(^l / / 

is monotonically decreasing and, provided it is bounded below, qualifies as a Lyapunov function 
guaranteeing the convergence of arbitrary initial data to the equilibrium solution. We note that the 
density g = is not constant since monomers can be added to or removed from the system. 

Finally there is a 'weak form' for the first equality in ([l|) 



r=2 



9lJl + l9r+l - gr]Jr 



(3) 



r=l 



In forthcoming analysis we assume initial conditions (cr(0)) which for large aggregation numbers 
decay faster than any exponential in r. 



II Coarse-graining of cluster aggregation space 

We now perform a coarse-graining contraction of the infinite set of Becker-Doring equations by 
systematically eliminating all the concentration variables except those which represent an aggregation 
number where 

A, = (r-l)A + l, r= 1,2,3,.... (4) 
2 



We then relabel the retained concentrations by The reduced fluxes are 

ar = TaA,aA,+i • • .aA,+i-i, (6) 

(3r+l = T6a,+1&A,.+2 • • • ^A.+i, (7) 

where T is a constant which represents a change of timescale; the kinetic equations then reduce to 

Xr = Lr-l — Lr, (r > 2). (8) 

This procedure is analogous to the Kadanoff block-spin renormalization procedure f^; detailed in- 
formation for cluster sizes between the aggregation numbers is lost. For more details of this 
procedure, see [|], |^. If the contracted system is to faithfully approximate the original system, we 
require that the special mathematical properties mentioned above are preserved under the coarse- 
grained rescaling. We can then draw on our renormalization procedure to extract the structurally 
stable phenomena present in the system. 

The physical properties of the full Becker-Doring system (|1|) are shared by the contracted system 
(11)^®- the partition function satisfies OrQ Ar = I^t+iQk+ii hence x^^ = Qa^Xi'' is formally an 
equilibrium solution. The function V{{xr}) = Jl"^! Xr{^og{xr /QAr^i'') — 1) has the same properties 
as (^. The weak form is still valid if Cr is replaced by Xr and Jr by L,.. Finally the density in 
the system is now defined by 

oo 

g = xi + Xj2[ir - l)X + l]xr. (9) 

r=l 

To apply renormalization ideas to this theory, we consider the repeated application of the coarse- 
graining transformation so we now reapply the contraction procedure with mesh size /x. 
Defining new variables zi = xi, Zr = and Ir as the flux from Zr to Zr+i, we find 

with Ar, Br determined from a^, (3r in an analogous way to (§|)-(|3)- A similar set of physical properties 
holds for this system of equations as for the original Becker-Doring equations. Thus a repetition of the 
coarse-grained contraction is identical to a single application with a larger mesh parameter A/i. This 
shows that it is sufficient to consider a system of equations which has undergone a single contraction 
with large A. 

Ill The case of constant coefficients 

Although ultimately a theory capable of handling arbitrary forms of rate coefficients and br+i is 
our goal, for the sake of simphcity let us start by considering constant coefficients - that is ar = a, 
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br = b. The parameter 6 = aci/b enables the system's behaviour to be classified. The cluster 
partition function is defined by Qr = {a/bY~^ and the forward coefficients in the reduced model by 
ar = Ta^, Pr+i = Tb^. Thus the size-independent rate coefficients ar = a, br = b are mapped to 
size-independent rate coefficients in the reduced model. This coarse-graining maps 6 to 6^, leading 
to three fixed points, ^ = 0, 1, oo. The large-time asymptotics of systems with constant coefficients 
have been analysed in detail in |^ , where it is shown that the 6 = case converges to the equilibrium 
solution, the 6 = oo case converges to the steady-state solution Xr = Xi by a diffusive wave which 
moves through aggregation space in such a way that its position is given by r = ~ t, and the 6 = 1 
case converges to the equilibrium solution Xr = Xi by purely diffusive means (x^ ~ Xierfc(r/2-\/t)). 

Ref 0] also shows that in order for the contracted system to preserve the correct large-time 
asymptotics, the parameter T in (|^)-(0) should take the value 

T= ^"^-^ - (11) 
A(aM-fe^)' 

this temporal rescaling implies that our renormalization is dynamic Following this temporal 
rescaling, the large-time limit of the density of the original system with 6 > 1 then scales with 
|ci(aci — which is identical to the result given by the coarse-grained system @; and both 
V{{cr}) of equation @ and \V{x) scale with — |ci(aci — &)^t^log6'. 

IV The case of power-law rate coefficients 

In many systems the reaction rates are not independent of size as assumed above, but rather depend 
on the size of the cluster according to some power law. We assume a,. = ar^, fe^+i = b^^, allowing 
us to model surface-limited aggregation in d dimensions with p = 1 — 1/d. The parameter 6 = aci/b 
remains a useful tool for classifying behaviour; the partition function remains Qr = {a/bY~^. The 
forward coefficients in the reduced model are 

ar = {[(r - 1)A + l][(r - 1)A + 2] . . . [rA]f . (12) 

For asymptotically large A these can be approximated by 

loga^ ~ Aloga + pA log(rA) — 1 + (1— r) log (^1 — ^) , (13) 

so for simplicity we shall take ar = (aX^r^)^, which is asymptotically correct at large r and differs 
only slightly at small values of r. In the same manner, the backward rate coefficients in the contracted 
model are given by Pr+i = {bX^r'^Y. 

Our coarse-graining contraction maps the set of models with power law rate coefficients into itself. 
The coarse-graining of power-law coefficients is only approximate when p ^ 0. However, the large- 
time asymptotics is qualitatively preserved, provided that a similar temporal rescaling is performed 
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as in equation (|rT|) 0. For any given model the contraction maps the exponent p to Xp. Following 
a contraction with large A, there are three cases to consider: p = 0, and large positive or negative 
p. The reduced system also has a different 6'-parameter, 9 = arX'^/Pr+i = 9^', thus the contraction 
maps 6 to 6^. The fixed points 6 = 0, 1, oo are therefore of most interest to us. Combining this 
information, there are nine fixed points of the coarse-grained contraction in {6,p) parameter space, 
and these form the basis of the ensuing analysis. 

Figure |I] shows schematically the effect of the contraction. In phase plane terminology, II has the 
form of an unstable node, I, III, V, VIII are saddle points (although they are at the limits of the 
allowable domain, so only have trajectories on one side of the fixed point), and IV, VI, VII, IX are 
stable nodes. Cases I, IV, VII all have partition function Qr = for r > 2; in cases II, V, VIII the 
partition function satisfies Qr = 1, whilst it is undefined in cases III, VI, IX since in all these cases 
the fragmentation rate is zero. Having no equilibrium configuration, these three cases approach a 
steady-state solution. 

V Effect of perturbations on the fixed points 

There are two reasons for wanting to study noisy coefficients: firstly any set of coefficients will be 
subject to uncertainties, whether derived from experimental data or a mathematical model. Secondly, 
systems are always susceptible to thermal (and in the models we study also spatial) fluctuations which 
locally alter the rate coefficients. In both cases it is necessary to know whether the models used are 
stable to minor variations in rate coefficients. 

Firstly we allow each reaction rate (a^, 6r+i for r = 1, 2, . . .) to be independently perturbed by a 
small amplitude random fluctuation of the form 

ar = arP{l + iy^r), K+i = br^il + i^Xr+i), r = l,2,.... (14) 

with z/ -C 1 and ^r, Xr+i being independent random variables with zero mean satisfying C,r, Xr+i = 
0{1). Such perturbations have no effect on the leading order equilibrium or steady-state solutions, 
or the large time asymptotics. 

A more interesting case is that in which the presence of noise in the rate coefficients is allowed 
to alter their leading order behaviour at large r. To examine these, we perturb the forward and 
backward rate coefficients according to 

Qr = ar^ + 6r, br+i = br^ + ^r+i, r = 1, 2, . . . . (15) 

where Sr,er+i have characteristic magnitude u <^1. We now investigate the effect of such perturba- 
tions on the equilibrium and steady-state solutions. Since we assume ci = 1, the partition function is 
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the equilibrium solution. However, as described in [Q, there are cases where the equilibrium config- 
uration formally has infinite mass and is hence not relevant; the system then approaches one of the 
family of steady-state solutions in which all fiuxes, J^, are equal. The steady-state fiux is determined 
by requiring the most rapid decay in as r — > oo. We now apply these ideas to the nine fixed points 
isolated earlier. 

Case I: p = 0, 9 = 0. Since the non-perturbed case has = for all r, the partition function and 
the equilibrium solution are then zero; introducing perturbations removes this degeneracy, and the 
equilibrium solution then becomes rapidly decaying in r, namely = Olu^'^^), where i> is a small 
parameter representing the typical size of perturbations 6r- 

Case II: p = 0, 6 = 1. In this case, introducing perturbations to the rates modifies the partition 
function from Qr = 1 to 

r 

1 + 5^(4- (16) 

k=l 

so small amplitude noise in the coefficients does not affect the leading-order behaviour of the system. 

Case III: p = 0, 6 = oo. In the absence of perturbations, there is no partition function for this 
case; when present, Qr ~ nfc=i(l/^A:)- However, this case converges to a steady-state rather than the 
equilibrium. When perturbations are included, the steady fiux is J = 1 + (61 — 62) + O^u^), which 
implies that the concentrations asymptote to = 1 + {Si—e2+er+i — Sr) + O^u"^). 

Case IV: p ^ 1, 6 = 0. As in case I, where noise is absent the partition function, Qr, is zero for 
r > 2. Introducing noise removes this degeneracy, for small u, Qr ~ 0{v'''^^). Thus, as in Case I, 
the equilibrium solution rapidly decays with r. 

Case V: p 1, 9 = 1. The balance of aggregation and fragmentation implies that = 1 in the 
case with no noise. The addition of noise to the rates alters this, to 

a = i + E(^^^)+^(^'). (17) 

wherein we see that the alteration to the partition function only affects the 0{v) correction term, 
leaving the leading order behaviour {Qr ~ 1) unaltered. 

Note that ii p > 1 then the system does not evolve to the equilibrium solution, but instead is 
attracted to a steady-state solution with more rapid decay in the limit r — 00. Perturbing the rate 
coefficients modifies this state to 

1 !^ 1 5k — Sk+i — Ji , . 

''^'-c(rt£^^? — — ■ <^'' 

where Ji = (1/C(p)) Hk'=ii.5k—^k+i)/k'^, which has constant fiux J = 1/C(p) + Ji- 
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Case VI: p ^ 1, 6 = oo. In this case the system approaches a steady-state solution, with flux 
J = 1 + (5i - 2-%2) + C(z/2), implying 

1 



1+ (5i-2-%2 + 7^^--) . (19) 
V (r + l)P rP)\ ^ ' 

Thus noise in the rate coefficients has a minor effect on the solution. 
Case VII: p <ti —I, 9 = 0. Formally we have, 

*^n^. (20) 

thus when r = 0{1), Qr = 0{u^^^). However, when r = Vc := 0{v^^^) 

Cr ~ [(z/i/P)!]-Pexp(z/i/nogz/). (21) 

For r > Tc, the perturbations have the same magnitude as the non-random part of the rate coefficient, 
thus all subsequent Qr values depend strongly on the perturbations 5fc,£^fc and have the order of 
magnitude given by (pOD- 

Case VIII: p <^ — 1, = 1. In the noiseless case this system converges to the equilibrium solution 
Cr = 1. When noisy coefficients are introduced, this solution may cease to be valid since at large r, 
the noise will be a leading-order effect. For small r we construct an asymptotic approximation to 
the modified equilibrium solution 

cr = l + E k-P{6k-ek+i) + 0{u^). (22) 
fc=i 

This approximation to the solution ceases to be valid at large r, where Cr = 0{1). We expect to 
remain 0{1) for all values of r, but to vary from c,. = 1 by significant amounts at large r. 

Case IX: p <^ —1, 9 = oo. In the absence of noise, this case approaches the divergent steady-state 
Cr = (with flux J = 1). For small amplitude noise, a modified form of this solution persists 



1+ 5i-2-%2- 



+ 0{u'y, (23) 



rP (r + l)P^ 

however, this ceases to be valid when r = 0{u^^^). For values of r of this magnitude and larger, 
perturbations to the rates cannot be neglected as they constitute a leading-order effect in the system; 
and Cr = 0{l/v) for all r > C(z/^/p). 

VI Effect of perturbations on the coarse-grained reaction rates 

In this section we examine the effect which the coarse-graining contraction procedure has on the 
perturbed rate coefficients. In particular, we investigate whether small amplitude noise in the full 
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description of the model maps to small amplitude noise in the reduced description. On inserting (^) 
into (D-(0) with T = A'^^, we obtain 

a, = aV^^ + A„ I3r+i = b^^^ + Er+u (24) 

where A^, Er+i represent the perturbations in the contracted descriptions and depend respectively 
on the Sk,ek- For each of the nine fixed points (in which a,b = 0, 1) we calculate the leading order 
form of this dependence. 

Case I: p = 0, 9 = 0. Since a = 0, 6 = 1, we have = A^. = O^u^), and Pr+i = 1 + O^u) Thus 
following contraction, the perturbations remain small. 

Case II: p = 0, 6 = 1. Following the coarse-graining contraction, the reaction rates are given 
by ar,(3r = 1 + 0{v). So the perturbations remain the same order of magnitude in the contracted 
model as in the full. 

Case III: p = 0, 6 = oo. The domination of aggregation is not altered by the presence of small 
noise, since ar = I + O^u) and Pr = Oiv). 

Case IV: p ^ 1, 6^ = 0. The contracted rates are given by ar = r^^ + Oiv) and I3r+i = 0{u)] in 
the latter, we have made the approximation ([T3|) valid for large r. The system remains fragmentation- 
dominated. 

Case V: p 1, 9 = 1. For large r, the rates in the contracted system have the form ar-, (3r+i = 
j-p>^ _l_ (9(z/) Thus the noise will not cause any change to the leading order form of the rate coefficients. 

Case VI: p ^ 1, 6 = (X). The domination of aggregation persists, since following contraction 
ar = + 0{v) whilst ^r+i = 0{u^). 

Cases VII-IX: p <^ —1. The formulas for Ar,Er+i in these cases are identical to cases IV,V,VI 
respectively. However, here p < so that at large cluster sizes r, the perturbations will be of the same 
order of magnitude as the deterministic part of the rate coefficients. This occurs when r = (9(z/^/p). 

In Cases I- VI, the noise indeed remains small in the contracted description of the model hence 
these may be termed universality classes, whilst in Cases VII-IX, this is not the case. In these last 
three cases, at large aggregation numbers, the noise in the full description is not small relative to 
the power law component of the rate coefficient, and this is reflected in the contracted model. In 
Cases VII-IX perturbations to the power-law rate coefficients play a major role in the kinetics at 
large cluster sizes r, as they do in the full model. Thus Cases VII-IX may be termed universality 
classes if the added noise decays faster than the given power law as r — >• cxd. 
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VII Conclusions 

We have applied renormalization ideas to the Becker-Doring model of cluster-formation. A novel 
feature of this work is that it is the cluster aggregation space which is rescaled, rather than a spatial 
dimension. Moreover, a dynamical renormalization is required to correctly maintain the time-scales 
of the growth and fragmentation processes following the rescaling of aggregation space. In the case 
of the power-law model, nine fixed points of the renormalization procedure have been identified and 
analysed in greater detail, quantitatively providing nine types of large-time asymptotics which may 
be exhibited by the system. Five of these systems tend to equilibrium, and the remaining four to 
steady-state solutions. The pure fragmentation cases (1,1V,V11) all tend to the trivial equilibrium 

In Cases 1, IV, Vll, a diffusive wavefront invades the large-r region where cluster concentrations are 
zero, leaving the equilibrium solution behind the wavefront. In Cases 11, V and Vlll the equilibrium 
solution is approached by purely diffusive mechanisms, no advection being present in the system. If 
p > 1 in Case V, then the system approaches a steady-state solution rather than the equilibrium 
solution, since the steady-state has faster decay at large aggregation numbers. This case is thus 
similar to Cases 111, VI, IX, all of which approach steady-states rather than true thermodynamic 
equilibrium. However, their large-time asymptotics are more akin to Cases 1,1V, VII, being dominated 
by a diffusive wave which moves into the large r-domain. 

Thus for the first time we have identified universality classes present in the Becker-Doring equa- 
tions, in that any system with power-law coefficients can be classified into one of the nine cases which 
correspond to fixed points of our contraction, and this qualitatively determines the system's large- 
time behaviour. In physical terms, our demonstration that a renormalization structure underpins the 
Becker-Doring equations carries with it the implication that universal behaviour can be identified 
in the approach of such systems to equilibrium or steady-states. In a forthcoming paper [Q, we 
shall discuss the temporal behaviour in detail. In the case of the Becker-Doring equations, this is a 
very welcome development, since it dispenses with the need to specify in full detail all the generally 
unknown fine-grained rate coefficients. For example the partition function is left unchanged by the 
coarse-graining, as is the equilibrium solution and the steady-state solution. At the end of Section 111 
we showed that the large-time behaviour of both the density and the Lyapunov function (free energy) 
were left invariant by our coarse-grain rescaling. It is perhaps worth pointing out here, however, that 
consideration of the asymptotic limit implied by the renormalization procedure is not necessarily 
always appropriate, e.g. for systems in which it is crucial to retain some level of fine-grained detail 
in order to properly capture the dynamics. 

The successful application of the renormalization techniques reported here opens the way for a 
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study of generalised Becker-Doring equations using similar methods; it also furnishes a firm theo- 
retical foundation for the analyses we have previously given of various generalisations of the basic 
Becker-Doring theory to a wide range of processes of physicochemical interest, including micelle and 
vesicle formation and self-reproduction 0, ^, generalised nucleation and growth phenomena 0, and 
macromolecular sequence selection in biopolymers []TU| . 
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